      function zlc = zpoly(n,z,nleja,zl)
% Calculates (z-zl_1)(z-zl_2)...(z-zl_n)  
 
      zlc= 1.;

      for ii=1:n
         zlc=zlc*(z-zl(ii));
      end

      end

